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^ ■ Abstract 



In 4D non-perturbative Regge calculus a positive value of the effective cosmological 
constant characterizes the collapsed phase of the system. If a local term of the form S' = 
'I2he{hi,h2,...} ^hVh is added to the gravitational action, where {hi, /12, ...} is a subset of the 
hinges and {A^} are positive constants, one expects that the volumes V/i^, V/ij, ••• tend to 
collapse and that the excitations of the lattice propagating through the hinges {/ii, /12, •••} 
are damped. We study the continuum analogue of this effect. The additional term S' 
may represent the coupling of the gravitational field to an external Bose condensate. 

04.20.-q Classical general relativity. 
04.60.-m Quantum gravity. 
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Since the first perturbative formulations of quantum gravity it was realized that the addition 
of a cosmological term / d'^x \J g{x) to the pure Einstein action gives the graviton a mass, 
which is positive if A < and negative - that means, the theory becomes unstable - if A > 

i- 

Nevertheless, whilst the vacuum fiuctuations of the quantum fields should in principle pro- 
duce a very large value of the cosmological constant, there is no observable hint of it. Namely, 
neither the Newtonian potential shows any finite range up to solar system distances, nor pure 
gravity exhibits any instability in the weak field case; finally the observations on cosmological 
scale set an upper limit on |A| as small as |A| < in natural units. This discrepancy 

is known as the "cosmological constant problem" and several possible escapes have been sug- 
gested 0]. Treating the problem at a fundamental level requires a non-perturbative approach, 
since one should be able to explain why the actual large-scale geometry of spacetime is fiat 
just from dynamic considerations. It is therefore not surprising that a definitive and generally 
accepted solution of the paradox of the cosmological constant is still remote. We shall not make 
any attempt to a new explanation here. 

We shall consider 4D pure gravity in the Euclidean approach, with special reference to quan- 
tum Regge calculus In this model the results of the numerical non-perturbative simulations 
compose the following picture of the behaviour of the cosmological constant: while its "bare" 
value A is generally nonzero, the effective value A depends on the energy scale ^ and vanishes 
at large distances like |A| ~ G'^^(/i/o)'^, where /q is the lattice spacing and 7 a critical exponent. 
This means that the quantum geometry fluctuates on small scales, but reproduces flat space at 
macroscopic distances. The sign of A is negative, thus the flat limit is well deflned. 

Our aim is to study the interaction of this gravitational system with a particular external 
source, namely a Bose condensate described by a scalar fleld (pi^x) = 4>o{x) + (t>{x). We assume 
that the vacuum density (f)Q{x) is forced from the outside to a certain value, as it can happen for 
instance in a superconductor subjected to external electromagnetic flelds. The positivity of the 
Euclidean action of ensures that the terms [d^(j)o{x)]*[d^(f>oix)] and m'^\(j)o{x)\^ act like positive 
cosmological contributions to the gravitational action, possibly inducing local gravitational 
instabilities. 

The structure of the paper is the following. In Section 1 we recall the main results of 
Euclidean 4D quantum Regge calculus concerning the effective cosmological constant. We 
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mention its scale behaviour and two possible interpretations of the lattice spacing Zq- In Section 
2 we consider in the continuum theory, at distances much larger than Planck scale, a weak 
Euclidean gravitational field g^v{x) = 5^i, + nhfj_^{x), whose action includes an infinitesimal 
effective cosmological term with A < 0. Such a theory can be regarded as the continuum limit 
of the lattice theory described in Section 1. We then couple minimally gfj,u{x) to the mentioned 
scalar field with vacuum density (f)o{x) determined from the outside, and write in detail 

the various terms of the action. It turns out that the infinitesimal graviton mass oc |A| 
receives a local negative contribution, which we denote by — This means, as we show 
explicitly for the simpler case of an almost-massless scalar field that if there exist some four- 
dimensional regions in which > m^, it is possible to find field configurations x{^) which 
make the Euclidean action unbounded from below. The field tends to develope singularities 
in those regions, or - if a cut-off mechanism comes into play - it tends to assume constant 
extremal values which are independent from those in the neighboring regions. We suggest that 
these "constraints" should be inserted in the equation for the propagator of the field. 

In Section 3 we remind the formula which gives the static potential energy in Euclidean 
quantum gravity as a functional average 0, ^. We illustrate its statistical meaning in the case 
of a weak field through an analogy with a simple 2D Ising model, which we also treat numerically 
with an elementary simulation. We show numerically that the insertion of local supplementary 
constraints due to an external field (in analogy to the singularities of the continuum case) 
damps in a sensible way the spin-spin correlations, and thus the interaction they represent in 
the model. 

1 Regge calculus. 

In the last years discretized quantum gravity on the Regge lattice has led, through the Mon- 
tecarlo numerical simulations of Hamber and Williams, to a better comprehension of the non- 
perturbative behaviour of the Euclidean gravitational field in four dimensions. Some of the 
features which emerge from their results depend on the specific model, but several oth- 
ers are quite general. In this approach the physical quantities are extracted from functional 
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averages. The partition function is written as 

d[g] e-'^"'^[^l (1) 



Geometries 

where the geometries are described by piecewise flat simphcial manifolds. The integral (i.e. the 
Montecarlo sampling) runs over the lengths of the links which define a field configuration. The 
action in (P has the form 

S = Jd'x^(^X-kR + ^aR^^p^R'^'P'^^ (2) 

or, in discretized version, 

_ r 4, A, 42^21 

(3) 



s[i]= E 

hinges h 



where Ah is the area of an hinge, Vh its volume and 5h the defect angle (see the original papers 
for the definitions of lattice quantities and the functional measure). 

In the following we shall set h = 1. The constants k and A are related, as "bare" quantities, 
to the Newton constant G and to the cosmological constant A: k corresponds to l/SvrG and A 
to A/SttG. It is important, however, to keep distinct the physical values G and A from k and A. 
The latter are entered as parameters at the beginning, and then a second order transition point 
for the statistical system described by Z is found by Montecarlo simulation. Actually, there is 
a line of transition, since one can also vary the adimensional parameter a, which does not have 
a macroscopic counterpart. On this line in the parameter space the theory admits a continuum 
limit. Unlike in perturbation theory, where a fiat background is introduced by hand, here the 
fiat space appears dynamically; namely, the average value of the curvature is found to vanish 
on the transition line, which separates a "smooth phase", with small negative curvature, from 
a "rough", unphysical phase, with large positive curvature. In this way the effective, large scale 
cosmological constant 



A = 4^ (4) 



vanishes in the continuum quantum theory. 

More precisely, the dependence of the effective (or "running") cosmological constant Ag// 
on the scale is the following. If we compute (^ on small volumes, the curvature fluctuates 
wildly. At larger distances the average curvature decreases, because the fluctuations tend to 



average out. If /x denotes the energy scale, close to the critical point the adimensional quantity 
|A|G behaves like 

(|A|G)(/i) ~ M\ (5) 

where Iq is the average spacing of the dynamical lattice and 7 ~ 1.56 in the first simulations. 
The sign of A is negative, as mentioned. Since the /^-dependence of G is quite weak, the effective 
cosmological constant decreases approximately like a power law as the length scale grows. 

Eq. (^) admits two different physical interpretations, depending on the role we attribute to 
the average lattice spacing Iq = ^ (P). In the usual lattice theories, /q is sent to zero in order 
to obtain the continuum limit. In this first interpretation, eq. (H) shows the way the lattice 
theory reproduces fiat space in the physical limit. No real physical meaning is assigned to the 
effective A which, at a fixed scale fi, is simply proportional to a positive power of the regulator. 

On the other hand, we can believe that in quantum gravity Iq has an intrinsic minimum 
value of the order of Planck length Lpianck ~ 10~^^ cm. This hypotesis arises independently 
from several operational models (for a review see ||^) or from more complex quantum theories 
(see for instance and references). In the framework of Regge calculus, it is possible to fix 
lo by imposing that the effective Newton constant G computed non-perturbatively is of the 
same magnitude order of the observed value; in this way one finds one more time that Iq has 
to be of the order of Planck length [|] . 

In this second case, that is Iq ~ Lpianck, the interpretation of eq. (|^) is different: it means 
that the effective cosmological constant tends to zero on large scale, while it is nonvanishing, 
in principle, on small scale (here and in the following we mean by "large" scale the laboratory 
or atomic scale, and by "small" scale the Planck scale). This interpretation does not neces- 
sarily have observable physical consequences, since in fact A could be far too small. Namely, 
the exponent gamma has been computed only for small lattices; an evaluation for lattices of 
"macroscopic" size is of course technically impossible, and thus only the experiments could tell 
us whether the law (^ keeps true for large distances, and with which exponent. The fact that 
on astronomical scale we have \A\G < 10~^^° for fi ~ 10~^^cm~^ constrains 7 to be approxi- 
mately larger than 2. But the vanishing could be much more rapid, so that we could disregard 
A at any physically relevant scale. 

Another property of Regge calculus which shows an intrinsic feature of quantum gravity is 
the instability of the "rough" phase with positive average curvature, i.e. with positive effective 
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cosmological constant. This phase does not admit any continuum hmit. Its fractal dimension 
is small, which denotes that the geometry is collapsed. 

The non-perturbative instability properties of the phase with positive average curvature 
extend the validity of known considerations based on perturbation theory. Namely, in the weak 
field approximation a positive cosmological constant in the gravitational lagrangian produces 
a negative mass for the graviton (compare and the next Section). Instabilities arise in the 
perturbative theory on a De Sitter background too [1^. On the contrary, a small negative 
cosmological constant generally does not imply any instability, but a small mass nig for the 
graviton, of the order of A^/^ in natural units. 

In conclusion, at a scale large enough with respect to the lattice scale it is possible to 
regard the discretized gravitational functional integral (|I]) as describing an almost flat mean 
field plus fluctuations. If the fixed point of the lattice theory is approached from the physical, 
smooth phase, the effective cosmological constant A (i.e. the average of the scalar curvature R) 
is very small and negative, and the large scale fluctuations of R are small too. The system is 
stable, because the field configurations with larger volume, in which the links are as stretched 
as possible, are preferred to the collapsed configurations, since the Euclidean action depends on 
the volume like 5* ~ / d'^x ^ ~ '^h- (On the contrary, a positive value of A would 

favour the collapsed configurations with smaller volume.) In this picture A can be regarded 
either as a purely formal regulator, which goes to zero in the physical limit /q 0, or as 
a physical quantity, though possibly extremely small (in the second interpretation, in which 

^0 ~ ^Planck)- 

Keeping these properties in mind, it is interesting to consider the case in which the coupling 
of the Euclidean gravitational field with an external source gives in some four-dimensional 
regions a positive contribution to the effective cosmological constant. In the next Section we 
shall analyse this phenomenon in the continuum case, that is, on large scale (in the meaning of 
"large scale" we precised above). 

On the lattice, such a coupling would correspond in the action (|^) to an additional term of 
the form 

E ^hVh (6) 

{hiM-} 

where {hi,h2...} is a subset of the hinges and Xh are fixed positive constants. 
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We expect that when the Montecarlo algorithm chooses for the random variation a hnk 
which ends in a hinge hie{hi, h2...}, the favoured variation will be that for which the volume 
Vh^ decreases. Thus the volumes {V/^^, Vh^, ■■■} will tend to collapse and the lattice excitations 
propagating through the hinges {hi, h2...} will be damped. 



2 Continuum case. 

We have seen that the results of quantum Regge calculus can be interpreted as leading at 
distances large compared to Lpianck to an effective Euclidean action for pure gravity of the 
form 

r . , r A 1 1 

(7) 



-R{x] 



_8nG SnG 

where the curvature fluctuations around fiat space are small and the effective cosmological 
constant A is negative and very small too. (As we saw in the previous Section, A is scale- 
dependent; we suppose here to stay at some fixed scale.) From the geometrical point of view, 
the small negative value of A stabilizes the system, preventing it from falling into small- volume, 
collapsed configurations. 

In the naive perturbation theory around the fiat background the A-term represents a small 



mass for the graviton 0. Namely, setting g^u{x) = 5^^ + i^h^i^^x), with k = v8nG, the determi- 
nant g of gfj^i, can be expanded as 

g = 1 + Kh^^^ + K^U^^ + ... , (8) 

where h^^\ h^'^\ ... denote expressions which are linear, quadratic etc. in h^^. The linear 
"tadpole" term h'^^'> is usually disregarded, since it is proportional to the trace /i^, which vanishes 
on physical states. The term h^'^'> takes the form of a graviton mass term, such that the mass 
is positive when A < (compare We thus have 



Seff = !d^x[ [mlU^\x) - i?(2)(x)] + A [/€/i(3)(x) + k'¥*\x 



+ 



[kR^^\x) + k'R^^\x) + ..]], (9) 



^As remarked in this widespread belief is not rigorously true. It has also been proved that all 
theories of a massive tensor field in Minkowski space which satisfy the usual QFT postulates are incompatible 
with General Relativity in the limit of vanishing mass. In our reasoning it is not essential to regard the 
cosmological term as a graviton mass term, but for simplicity we stick to this terminology. 
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where nig is proportional to A^/^ and R^'^\ R^^^ ... denote the parts of R which are quadratic, 
cubic etc. in h^^. 

The first bracket contains the quadratic part of the action. The third bracket contains the 
famihar self-interaction vertices of the graviton, involving respectively 1, 2, ... derivatives of 
h^i^. The second bracket contains self- interaction vertices which are peculiar of the theory with 
A 7^ and do not involve derivatives. We are however not interested in the self-interaction 
vertices of hf^i, in the following. 

Now we would like to consider the interaction of h with a scalar field having non- 
vanishing vacuum expectation value 0o- We suppose 0o to be spacetime dependent and denote 
(f){x) = (f)o{x) + 4>{x)] (f)o{x) is regarded as a quantity determined from the outside, that is, 
as a source term. In this way, the field (f){x) describes a Bose condensate with ground state 
density (poi^) fixed by external conditions (like, for instance, the Cooper pairs density in certain 
superconductors under external e.m. field). 

The total action is 



S = Seff + J d^x [d^ [00 (x) + 0(x)] * [<Po{x) + ^{x)] g^^x) + ^m^ |0o(x)| 



(p*o{x)4>{x) + 0o(x)0*(a;)] + ^m^\^{x)\^} 



(10) 



We can rewrite it as 



where 



S 



d^x \[g{x) 



^ 1 2/ A 



-^i?(.)} + ^. + ^. 



'11^ 



^^\^) = -[9^0S(^)][5>o(x)] + -mJ|0o(x)r; 



^1 

^2 



2 

- / d'^x^Jg{x)^^,(t)l{x)^^(t)Q{x)Kh^'''{x) 



(12) 
(13) 



- j d^x ^ g{x)[ml\4){x)\^ + ml [0S(x)0(a;) + 0o(a;)0*(a;)] + 

+ [d,4>*{x)d,~^{x) + d,cl>l{x)d,~^{x) + d^4>*{x)d,Mx)\ g'^'^ix)} (14) 



For brevity, we have not expanded here ^ and R like in eq. (P). Let us first look at 
the terms Si and 5*2. The term Si describes a process in which gravitons are produced by 
the source ip^ix). The term 5*2 contains the free action of the field 0(x), which describes the 
excitations of the condensate, and several vertices in which the graviton field h^y{x) and 0(x) 



interact between themselves and possibly with the source. All these interactions are not of 
special interest here and are generally very weak, due to the smallness of the coupling n. The 
relevant point is that the purely gravitational cosmological term receives a " local" positive 
contribution |yU^(a;) which depends on the fixed external source (f)o{x). According to eq.s (y), 
(|^ and to our discussion of the sign of A, this amounts to a negative mass contribution and 
could lead to instabilities. 

Let us study the effect of such a local negative mass contribution in the simpler case of a 
scalar field in fiat space. We consider a scalar field x with very small mass and add to its 
free Euclidean action a source term of the form — |/i^(a;)x^(a;), which represents a "localized 
negative mass" . The action becomes 

Sx = lJ d'x{[d^x{x)][d,x{^)]+mlx\x)-f,\x)x\x)} (15) 

To fix the ideas, let us suppose that is different from zero in certain four-dimensional 

regions Qi, where it takes the constant values /ij. Outside these regions let go rapidly to 

zero. The solution of the classical field equation for x is obtained by minimizing the action. 
In the absence of the source term we would have of course x{x) = const. = 0, because this 
minimizes at the same time the gradient term ^[df^x{x)dfj_x{x)] and the mass term im^x^(x) 
(both positive-defined). 

In the presence of the source term — i/i^(x)x^(a:), the action decreases when becomes 
large within the regions Qi. The growth of x iii these regions is limited only by the fact that 
outside them x must go to zero (due to the term m^) and that the gradient in the transition 
region cannot in turn be too large. Let us suppose for instance that there is only one region 
Qq, with the shape of a 4-sphere of radius tq. Let /i(x) take the value /xq inside Qq and zero 
outside. We consider for x the trial function x(r) = Xo /('")• The action becomes (we admit 
spherical symmetry) 

= ^'Xl f drr' {[/'(r)]^ + mj[/(r)]^ - /.^^(r - ro)[/(r)]^} (16) 

We see that if the integral is positive, the value of xo which minimizes the action is still xo = 0. 
On the contrary, if the integral is negative, the action is not bounded from below as xo grows. 
We choose the following explicit form of /(r): for r < Tq let /(r) = 1, i.e., x(^) = Xo] for 
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r > tq let be /(r) = exp[— (r — rQ)/5]. We thus have 

= vr^X^{(^ + mj)£°°cirr^e-(-o)/^-l(/.g-mjK^} (17) 
= { (1 + m^) (^^^ + ^ro5^ + \rl5^ + '-r^s) - ^(^^ - mJ)ro^} (18) 

It is easy to check that for suitable choices of the parameters the expression within the braces 
in eq. ([TsD turns out to be negative. For instance, if can be disregarded with respect to /io 
and 6 ~ tq, the expression is negative provided /iq ^ S^^. If 5 ^ tq the expression is negative 
provided yUo ^ rQ^, etc. Thus the system is unstable. For the graviton the instability is even 
worse, because the kinetic term R^"^^ is not positive-definite. 

Physically, we might of course invoke some additional "regularizing" process which comes 
into play for large values of xo and cuts the action. Thus the value of x^{^) inside the region 
is forced by the source to a certain maximum, and this will affect the propagation of the field. 
The differential equation for the propagator of xi^) in the presence of the source must now 
satisfy additional boundary conditions on each region fij. This amounts to a very complex 
mathematical problem; the physical consequence will be a "damping" of the correlations of 
x(x). One can check this numerically in some model; we shall do this shortly in the next 
Section, referring to a toy bidimensional spin model. 

3 The formula for the static potential energy. 

In this Section we remind the formula which gives the static potential energy in Euclidean 
quantum gravity as a functional average 0, ^. We illustrate its statistical meaning in the case 
of a weak field through an analogy with a simple 2D Ising model, which we also treat numerically 
with an elementary simulation. We show numerically that the insertion of local supplementary 
constraints due to an external field (in analogy to the singularities of the continuum case) 
damps in a sensible way the spin-spin correlations, and thus the interaction they represent in 
the model. 

Let us first consider, in Euclidean quantum field theory, a scalar field $(x) with action 
5'o[$] = / (i^x L($(x)). In the presence of the external source J{x), the ground state energy of 
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the system can be expressed as 

n Jdmexp{~h-^ [J d^x imx)) + J d^x^{x)J{x)]} 

E = lim ——log 5^ 7 ^ 

T Jd[^]exp[-h-^ Jd^xL{^{x))} 

= lim ~\og(expl^-h-^ I d^x^ix)J{x)^y (19) 

where it is assumed that the source vanishes outside the interval (— T/2,T/2) and that the 
couphng between J and $ is hnear. More generally, a formula similar to (p]9| ) holds when we 
are dealing with more fields $a and corresponding sources J^, and when the coupling between 
fields and sources is not linear. 

As an useful application of eq. ([T^, we can write the interaction energy E{L) of two static 
pointlike sources of the field, kept at a fixed distance L. We just need to insert the suitable 
expression for J. The trajectories of the two sources are in fiat space 

x^.ih) = (ti, 0, 0, 0); x^(t2) = (t2, L, 0, 0). (20) 

In ordinary gauge theories we may re-obtain in this way the Wilson formula for the static 
quark-antiquark potential. In quantum gravity we are led to the following equation for the 
static potential 

^ J d[g] exp l^-h'^ S[g] + Ei=i,2miJ\ dt ^ g^^\xi{t)\x'^{t)x'({t) | 
= T^™ -T^°^ ^d[g] exp{-n-'S[g]] 

(21) 

= lim — — log ( exp < —h^^ I ' dsi) ) (22) 




S 



where S is the euchdean action 

In eq.s (^l]), ( ^2|) the lines Xi{t) and X2{t) must be parallel with respect to the dynamic metric 
gfj_u and thus they should in principle be re-traced for each field configuration of the functional 
integral. In practice, it is extremely difficult to compute a functional integral defined in such a 
formal way. Let us then limit ourselves to consider weak fiuctuations of the gravitational field 
about fiat space. The trajectories of the sources may be defined with respect to the background 
^Notice that this formula apphes also to two masses which are not pointlike ("pointlike particle" is actually 
an ill-defined concept in General Relativity), provided we can disregard the internal degrees of freedom. Namely, 
their action is still equal to J2i I dsi, where the integrals are taken along the trajectories of the centers of mass. 
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metric like in (pOl). It is straightforward to reproduce in this way the Newton potential energy 
Q ; higher order corrections have been computed too ||T2| . The geodesic distance between the 
trajectories xi{t) and X2{t) is now equal to L only on the average; in fact such an approximation 
is not without physical meaning, since in any realistic source the fixed distance at which the 
two masses are kept can only be an average value. Also in the non-perturbative evaluations 
of eq. (p2D in quantum Regge gravity p , the distance L is evaluated a posteriori as the mean 
value of the geodesic distance on all configurations. 

Let us consider the almost-flat metric g^^, = + rj^y in a fixed gauge. We obtain 

rT/2 rT/2 . 1 i-T/2 

/ rfs, -/ dtiJl + h^^[x{U)]c^T + - dUh^^[x{ti)] (23) 

J-T/2 J-T/2 ^ Z J-T/2 

and we see that to a first approximation the fluctuations of h have the effect of making each 
line "shorter" or "longer". Let us call a[h] the (gauge-invariant) difference between the lenghts 
of the two lines in a field configuration h, and assume for simplicity that the masses of the two 
sources are equal: mi = m2 = m. We may expand the exponential in eq. (P^, finding (note 
that obviously vanishes by symmetry) 

E = 2m+ Yirn^-- log (e^v{-h-^ma[g\]) = 2m + ^lim (a'lg]) + ... (24) 

This equation exibits an interesting relation between the vacuum fluctuations of the geometry 
and the static gravitational potential. To illustrate better its "statistical" meaning we would 
like now to introduce a toy analogy with the 2D Ising model. 

Let us consider a planar spin system with periodic boundary conditions and the local cou- 
pling H = —JJ2i,j SijSi'jf (J > 0; neighbours of Let us then consider two columns 
ji and j2 at a distance of L lattice spacings (see fig. 1). 

We can regard this system as the analogue of a discretized configuration of a 4D gravitational 
field on the plane between the two parallel lines of eq. (pOl). The spin variables ±1 represent 
fluctuations of the metric. At the transition temperature, the fluctuations of the spin variables 
along the two lines are correlated, approximately like 1/L. 

Going back to eq. (p3| ) and making the correspondence s ^ hu, we see that the analogue 
of a[g] is the difference between spin sums taken along the columns ji and j2- The analogue of 
eq. dl) is 

12 



Siji 1 + ( ^ Sij^ \ 2 ^ Sij^ ^ Sfcj2 
i / \ i / i k 



(25) 



We are interested only in the term which depends on the distance L between the two 
columns, that is, the product term 

Pi2 = (l^Siiil^Sfcj^)- (26) 
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Fig. 1 - Spin sums taken along two columns of a 2D Ising system. 

A numerical simulation with a simple 10 x 12 system has given, as expected, the following 
results at the critical temperature: 

L (lattice spacings) p_{12} 



2 


1. 


.15 


+/- 


0. 


,13 


3 


0. 


,74 


+/- 


0. 


,09 


4 


0. 


,48 


+/- 


0. 


,09 


5 


0, 


,42 


+/- 


0. 


,08 


6 


0, 


,41 


+/- 


0, 


,08 



We see that the L ^ law is approximately verified also for the correlation between the spin 
sums taken along the two columns. Only for L = 6 there is a deviation, which can be explained 
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as due to the periodic boundary conditions. 

We observe that, due to the nature of the spin system, the correlation between two spin 
sums along ji and j2 is necessarily positive. Thus, the part of {cijsing) which depends on the 
distance L is negative. In perturbative quantum gravity one finds instead that the correlation 
{hii{x)hii{y)) is always negative, which leads to the correct negative sign for the potential 
energy in eq. (|^). 
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Fig. 2 - Insertions of singular points "*" between the columns 
jl and j2. 

We now introduce some supplementary conditions, in order to simulate the case in which the 
spin variables assume on certain sites a fixed value. This could be due, like in the gravitational 
case to which we are interested, to the localized action of an external field. In the spin model 
we may imagine that an external magnetic field localized on certain sites forces spin-flips. With 
reference to Fig. 2, let us suppose that the spin-flips occurr on the sites marked with a star 
and placed between the two columns ji and j2 whose correlation we are measuring. (The two 
columns are denoted by "x", while all the remaining sites are denoted by "0".) To prevent an 
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uncontrolled "driving" of the total magnetization, we associate to each spin- flip an opposite 
flip on a neighboring site. The flips occurr at each Montecarlo step; since the mean frequency 
at which the regular sites are flipped during the simulation is 120 times smaller, the resulting 
effect is to force a zero at the "*" sites. 

The precise positions of the sites at which the spin-flips happen are almost irrelevant; we 
find in all cases, as appears from the following table, a sensible diminution of the correlations 
between the spin sums taken along the columns ji and j^: 

Number of "*" sites p_{12} for L=5 



3 (Fig. 2. a) 0.19 +/- 0.04 

4 (Fig. 2.b) 0.13 +/- 0.03 



It appears therefore that the insertion of variables which are driven by an external field 
damps the correlations in the system, and that this mechanism is of a quite general nature, 
although we are not able to give a precise analytical description yet. 

4 Conclusive remarks. 

We have investigated in this paper an unusual interaction mechanism between gravity and a 
macroscopic quantum system driven by external fields. This idea was originally suggested by a 
possible phenomenological application O, but the mechanism is interesting also from the 



purely theoretical point of view and deserves further numerical and analytical investigation. 

We have shown that under certain conditions the gravitational field becomes unstable and 
may develope singularities, but we have not tried to find a physical regularization and to 
compute the effect of the regularized singularities yet. Simple physical analogies show however 
that they generally reduce the gravitational long-range correlations. Our next task will be the 
esplicit estimation of the changes in the correlation functions in terms of the squared density 
|0o(a;)P of the Bose condensate and of its squared gradient [df^(po{x)]*[d^(po{x)]. 
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